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Abstract 

We study the properties of superconductor in presence of a finite concentration of repulsive centers. The supercon- 
ductor is described by the negative U Hubbard model while repulsive centers are treated as randomly distributed 
impurities with repulsive interaction. Analyzing the paring potential and local density of states at impurity sites 
we find a wide range of the system parameters where the it - like state could possibly be realized. Comparison of 
our results to the single repulsive center case is also given. 
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A superconducting system is regarded as in the ir- 
phase if there is a sign change of the order parame- 
ter between two subsystems. The simplest example is 
the junction made from two superconductors with the 
phase of the order parameter equal to 7r [1,2]. In this 
case Josephson current becomes negative in contrast to 
the usual 0-phase junction. Another example are gran- 
ular high-T c materials which can likely form network 
of microscopic n -junctions [3] between small regions 
with different phases of the order parameter. In such 
systems the zero-energy Andreev bound states, zero- 
bias conductance peaks, paramagnetic Meisner effect 
and spontaneously generated currents take place [1,2]. 
In our recent report [4] we have considered a single 
impurity with repulsive interaction embedded in a s- 
wave superconductor on two dimensional lattice show- 
ing conditions which should be fulfilled to obtain the 
'it' state. For a uniform embedding system such a state 
has been found in the limit of very strong repulsion at 
the impurity site. In the present paper we discuss an- 
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other possibility. Our superconductor possesses a finite 
concentration of repulsive centers with random distri- 
bution. With such assumptions V sates could appear 
for moderate values of repulsive interactions at impu- 
rities due to additional disorder in paring potential A;. 

The system is described by random U Hubbard 
model [5] with the Hamiltonian 



+ 



(1) 



where i, j label sites on a square lattice, tij = —t 
is the hopping integral between nearest neighbor sites 
and n is the chemical potential. Disorder is introduced 
through the random site interaction Ui. Depending 
on site i (i = A or B with binary alloy distribution 
Ai_ c B c ) Ui describes attraction (Ua < 0) or repulsion 
(Ub > 0) between electrons with opposite spins occu- 
pying the same site. In the Hartree-Fock approxima- 
tion [6] (we dropped the Hartree terms Ui{rii a ), and 
considered only the half field system n — 1 , for simplic- 
ity) the corresponding Gorkov equation at zero tem- 
perature (T — 0) has the following form 
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Fig. 1. Superconducting order parameter A versus concen- 
tration of repulsion impurities c for Ua/W — 0.3 (a) and 
Ua/W — 0.3 (b). Local density of states Nb(Ef) at the Fermi 
energy (c) for U A /W = 0.6 and U B /W = 0.3, 0.6, 1.0. 
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where 



A, = UiXi = —Ui— I duImG 12 (i,i;ui+i0), (3) 

7T ' 



and Ef = fJ.(T = 0) is the Fermi energy. 

Using the Coherent Potential Approximation (CPA) 
to treat disorder in the paring potential Aj [5] we have 
found a wide region of system parameters c, Ua, Ub 
where V state could exist. In that cases Aa and As 
were of different sign and Nb (Ef) 7^ 0. In Fig. la and b 
we present the superconducting order parameter A = 
(1 — c)Aa + cAb for attractive interaction Ua = 0.3 
(Fig. la) and Ub = 0.6 and a number of Ub repul- 
sions. We have checked that in the first case Ua = 0.3 
we have clear gap and Nb(Ef) — up to some critical 
concentration of c = Co, above which the system be- 
comes normal A = 0. In the second case Ua = 0.6 the 
systems stays superconducting for any c but transits to 
'it' state for a critical value of concentration c = ci de- 
pendent on Ub ■ Larger Ub leads to smaller c\ . Above 
this concentration Nb(Ef) ^ (Fig. lc). 

For better clarity in Fig. 2a,b we have plotted local 
density of sates for Ub = 0.6 and Ua = 0.3 (Fig. 2a) 
and Ua = 0.6 (Fig. 2b) Note that in both considered 
cases (Ua/W = 0.3 or 0.6) the impurity bands are 
symmetrically located with respect to the Fermi energy 
as it is in case of a single repulsive impurity [4]. 
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Fig. 2. Local density of states (LDOS) at the B site for 
U B /W = 0.6 and U A /W = 0.3 (a) or U A /W = 0.6 (b). 



Our results in Fig. 2a show transition from super- 
conducting system c = 0.1 and c = 0.3 with a clear gap 
to normal state c = 0.7. This confirms the results with 
a single impurity [4] where it has been shown that for a 
relatively small bulk interaction the 'it' state is absent. 
On the other hand Fig. 2b shows the evolution of sys- 
tem with a clear gap around Ef (up to c = 0.3) where 
N B (E) yi for any E but not E — Ep (N B (E F ) = 0). 
Finally for c = 0.7 Nb(Ef) / which could possi- 
bly correspond to a state with 7r-like properties. Ac- 
count for relatively large values of the random poten- 
tial Ui make Hartee-Fock approximation questionable 
so the results obtained here should be treated as qual- 
itative. The presence of Hartree terms Ui{nicr) could 
also change the results in some range of impurity con- 
centration c forcing the transition to normal state for 
a large enough interaction Ua [5]. 

In summary we have studied properties of the su- 
perconductor with repulsive centers and showed that 
for a large concentration of those centers and a large 
on-center Coulomb interaction system could possibly 
evolve to the n state. Two necessary conditions for ex- 
istence of the 7r state, namely, the sign change of the 
pairing amplitude and non-zero the density of states 
at the Fermi level are fulfilled. 
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